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Bubbles display astonishing acoustical properties since they are able to absorb and scatter large 
amounts of energy coming from waves whose wavelengths are two orders of magnitude larger than 
the bubble size. Thus, as the interaction distance between bubbles is much larger than the bubble 
size, clouds of bubbles exhibit collective oscillations which can scatter acoustic waves three orders 
magnitude larger than the bubble size. Here we propose bubble based systems which resonate at 
frequencies that match the time scale relevant for seismogenic tsunami wave generation and may 
mitigate the devastating effects of tsunami waves. Based on a linear approximation, our naive 
proposal may open new research paths towards the mitigation of tsunami waves generation. 



Tsunami waves are one of the most devastating natural 
events. Recent examples originated by submarine earth- 
quakes having a magnitude Mw > 9.0 have been cov- 
ered by worldwide media showing the outcomes of such 
a natural disaster and the human and economic tragedy 
thereafter. Current measures and efforts are focused on 
forecast, early warning systems, and occasionally also 
on coastal tsunami run-up mitigation [1]. Paradoxically, 
such destruction is generated by a very small amount 
of the strain energy released by the faulting, roughly 
speaking less than 1% [2, 3]. Then, one may ask: Could 
it be possible to mitigate the generation of seismogenic 
tsunami waves?. Any thinkable mechanism that can be 
thought for this purpose acquires enormous proportions. 
Mitigation of coastal waves would imply huge barriers 
over extensive zones that are not affordable in the case 
of large tsunamis. It has been proposed the deployment 
solid periodic resonators arrays as a way to block water 
waves [4]. However, if applied to tsunami scales it con- 
cerns very large rigid scatterers hardly feasible for real 
situations. 

In a similar way to what happens in the sky, where kilo- 
meter size clouds composed by water droplets can be seen 
thank to light being scattered by tiny water droplets, we 
propose the use of air bubble metaclouds immersed into 
water to affect the generation of tsunamis (see Fig. 1). 
Thus, instead of dealing with a huge system composed 
of large parts, our proposal takes advantage of properly 
arranged small building blocks that could collectively res- 
onate at a frequency scale a tsunami requires. 

Bubbles have unique properties due to the large mis- 
match between air and water mechanical compressibil- 
ities [5, 6]. A bubble (Fig. 1 stage 0) having equilib- 
rium radius Ro displays a resonance [7] whose angular 
frequency is given by ujq = Rq 1 y/S'ypo/ where p corre- 
sponds to water density, 7 represents the specific heats 
ratio for air, and po denotes the static pressure in wa- 
ter. Considering the sound speed in water as c = 1500 
m/s and atmospheric pressure, one can easily find that 
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Figure 1. Diagram depicting the metacloud approach. 

The building block (stage 0) consists of small air (gas) bubbles 
of radius Ro which arranged form a spherical cloud (stage 1) 
of radius R\ with a gas filling fraction a±. A number of iV2 
clouds are then arranged within a disk (stage 2) of radius R2 
and thickness 2Ri with a cloud filling fraction forming a 
metacloud. From the resonance formulae, n, j = 1,2, . . . . 



\o/2Ro « 200, which means that the wavelength in wa- 
ter at the bubble resonance is 200 times larger than its 
size. This deep subwavelength behavior makes bubbles a 
perfect building block for a metamaterial [8] . In order to 
figure out the length scale of a bubble resonating with a 
tsunami wave, we should estimate the order of magnitude 
of the time scale of the seismic movement responsible of 
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the tsunami formation. Long gravitational surface waves 
(tsunami waves in open ocean) in the linear regime prop- 
agate with a speed c g = ^fgh over the water layer of 
depth h (g = 9.8 m/s 2 is the acceleration due to grav- 
ity). Choosing the time interval that this wave requires 
to propagate a distance equal to h as a target we obtain 
frequencies on the order of uj^/2tt = ^g/h = 0.075 Hz 
for h = 5 km. Using the resonance frequency of an 
hypothetical bubble tuned at ujt we obtain a radius 
i?T = ^ a/ S'JPo I p ~ 43 m. Although having a deep 
subwavelength resonance, this enormous bubble suffers 
from several practical drawbacks such as buoyancy that 
may be prevented by using huge gas balloons attached to 
the see bottom. However, the gradient pressure appear- 
ing at increasing sea depths should strongly affect the gas 
balloon sphericity, modifying its scattering properties. 

Figure 1 shows the strategy we have followed to reach 
a resonance frequency of 0.075 Hz using a small bubble of 
few millimeters as building block (stage 0). In the next 
stage (stage 1) a bubble cloud is considered and even- 
tually the metacloud stage (stage 2) is reached. Bubble 
clouds can be found in nature and have been studied as 
a source of noise in underwater environments [9, 10]. In 
this context, collective modes of bubbles in clouds (Fig. 1 
stage 1) [6] have been observed due to the large acoustical 
scattering and absorption cross sections (<j s , <j a respec- 
tively) of single bubbles allowing long range interactions 
across the cloud. 
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Figure 2. Comparison between the interaction radii of 
a cloud and bubbles of different sizes. Interaction radius 
^ct/tt in meters calculated from (a) scattering a s and (b) 
absorption a a cross sections for a cloud of radius R\ — 37 m, 
a large bubble of the same radius Ri, a bubble of radius R g 
having the same volume of gas than the cloud, and a bubble of 
radius Rq as a function of Ri/X. The cloud has Ri/Ro — 10 4 
and oli = 0.024. 

Defining the acoustical interaction radius as \J 'a /ir, we 
can compare between a spherical cloud and bubbles (or 
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Figure 3. Variation of the interaction radii as a func- 
tion of the cloud radius. Maximum acoustical interaction 
radius for (a) scattering and (b) absorption as a function of 
the cloud radius Ri and compared to that of a bubble of the 
same size R±, a bubble of radius R g having the same gas 
volume than the cloud, and a bubble of radius Ro, so that 
R 1 /R = 10 4 and ai = 0.024. The vertical arrows indicate 
Ri = 37 m, for which the interaction radius are plotted in 
Fig. 2. 



balloons) of different sizes as in Fig. 2. Following [11] and 
considering thermal effects and damping mechanisms, we 
calculate scattering and absorption interaction radius in 
the linear regime for time harmonic excitation, i.e. the 
bubble radius can be written as R = Rq(1 + (pexp(iujt)) 
when the complex oscillation amplitude \ip\ « 1. Then, 
the bubbly water of the cloud [12], can be treated as an 
effective homogeneous medium coupled to the surround- 
ing liquid. A cloud of Ri = 37 m formed by bubbles 
of Ro = 3.7 mm with a gas filling fraction a\ = 0.024 
clearly shows a collective resonant mode at R\/\ ~ 0.01 
(Fig. 2 see labels). As it can be expected, the cloud inter- 
action radius curves are orders of magnitude away from 
that of their building block. Also, the maximum inter- 
action radius of the cloud overcome those provided by 
large bubbles having either the same radius or a ra- 
dius R g = 10.6 m which contains the same volume of gas 
than the cloud. Although both large bubbles have lower 
resonant frequencies, the interaction radius is larger for 
the cloud, which depends strongly on the radius R\. Tak- 
ing (y^cr/^Max for different cloud radii R\ while keeping 
Ri/Ro = 10 4 and the same gas filling fraction for the 
cloud (ql\ = 0.024), we have plotted in Fig. 3 the max- 
imum acoustical interaction radius of a cloud as a func- 
tion of R\. For the sake of comparison and in the same 
manner as in Fig. 2, we also depict the cases of a single 
bubble (or balloon) of different sizes. It can be seen that 
above a critical radius {R\ w 10 m) the cloud has higher 
interaction radius than the corresponding large bubbles. 
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Figure 4. Metacloud resonances. Real part of the time harmonic component of the clouds pressure Re{0 n } (color scale and 
cloud size) in a metacloud for (a), (d) square and (b), (e) random metaclouds of radius R2 = 1.7 km and a cloud filling fraction 
a\ — 0.3 and N\ — 1005. (c), (f) \<j> n \ in normalized units as a function of r / R2 for both square and random arrays. The top 
row panels correspond to the first mode at uj\/2ti = 0.075 Hz and the bottom row ones to the second mode at 002 /2tv = 0.15 Hz. 
The clouds sizes are exaggerated to allow a proper visualization. 



The acoustic results shown above allow us to consider 
low frequency modes of clouds of bubbles as a single bub- 
ble, i.e. the lowest resonant mode of a cloud is similar to 
the fundamental breathing mode of a single bubble. Also 
the collective mode takes place when the wavelength in 
water is at least 50 times larger than the cloud radius 
Ri — 37 m [12], which corresponds to a length around 
1.85 km. This wavelength value can further be enlarged 
to the tsunamis length scale by increasing the size of the 
bubble cloud. Alternatively, we can move further up to 
larger scales using the spherical cloud as a building block 
for a two dimensional (2D) metacloud (stage 2 in Fig. 1) 
and expect to see similar effects. In fact, using the multi- 
ple scattering method [13] [14] we are able to observe the 
interaction between clouds within a metacloud, as de- 
picted in Fig. 4. Expressing the pressure at the n-cloud 
as p n = po + (j) n exp(zcjt), then Re{(j) n } will give us infor- 
mation concerning the pressure at t = as depicted in 
Fig. 4(a), (b), (d), and (e). Fig. 4(c) and (f) show the 
time averaged pressure at the n-cloud \<fr n \ as a function 
of its distance to the metacloud center normalized by the 
metacloud radius r/i?2- These two modes of the meta- 
cloud are obtained arranging N± = 1005 clouds to form 
either square or random two-dimensional arrays of radius 
R2 = 1.7 km having a cloud filling fraction a 2 — 0.3. 
Both modes have the same shape regardless the way the 
clouds are arranged (Fig. 4 (a) and (b); (d) and (e)), 



which can be directly observed in Fig. 4 (c) and (f). In 
addition, there is a good qualitative agreement between 
the mode shape for spherical bubble clouds, for which 
(j)(r) = sin(kr)/kr as derived in [15] (k is the wavenum- 
ber in the bubbly effective medium), and present results 
for a circular metacloud as shown in Fig. 4(c) and (f). 
Thus, our results predict the existence of metacloud col- 
lective modes starting from large bubble clouds whose 
acoustical characteristics are comparable or even better 
than its large bubble counterparts. Consequently, a suit- 
able spatial distribution of tiny bubbles having 3.7 mm 
in size would be able to collectively oscillate at the low 
frequencies (0.075 Hz) of tsunamis. 

In summary, within the linear approximation and us- 
ing the acoustical approach, we predict the existence of 
collective oscillations of metaclouds. This phenomenon 
could be used to scale down the resonant frequency of a 
system able to target time scales which are characteristic 
in the generation of seismogenic tsunamis. Our proposal 
concerns a naive approach as it avoids important ingre- 
dients such as bubble cavitation and nonlinear effects in 
clouds. An extended discussion on the role of incompress- 
ibility in the formation of seismogenic tsunami waves and 
metaclouds, the nonlinear behavior of bubbles and bub- 
ble clouds as well as some notes on the feasibility of the 
metacloud approach are given in the Appendix. We hope 
this work could open a new path towards tsunami gener- 



4 



ation mitigation, which, to the best of our knowledge, is 
not yet included in any tsunami-related agenda[l]. Fur- 
thermore, our study could stimulate further investiga- 
tion towards experimental demonstration of the collective 
oscillations of bubble metaclouds and the development 
of more realistic and quantitatively accurate theoretical 
models that may lead to a feasible strategy to mitigate 
tsunami generation. 
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Figure 5. Power balance of the two-dimensional metacloud 
(i?2 = 1.7 km, Ri — 37 m, and Ro — 3.7 mm) calculated 
using multiple scattering method under Born approximation. 



Appendix 

Born approximation The effects of the bubble cloud 
on the acoustic wave propagation are difficult to deal 
with even within the linear regime. Born approximation, 
which is assumed in the multiple scattering formulation, 
limits the amplitude of the scattered wave to be negligi- 
ble in comparison to the incident wave. In bubble clouds, 
however, this is not always the case and the model fails 
in giving quantitative predictions even at low filling frac- 
tions (a\ ~ 0.01) [16]. The power balance of the sum of 
the scattered and the absorbed power over the incident 
power IFr/IIi deviates from unity as the Born approxi- 
mation fails. However, the position of the peaks where 
the power balance is not preserved corresponds to the fre- 
quencies where collective modes appear. We have tested 
this behavior for clouds and expect the same to hold for 
metaclouds (see Fig. 5). Thus, we obtain a qualitative es- 
timation of metacloud collective modes depicted in Fig. 4 
of the paper. More sophisticated methods have been de- 
veloped for the one-dimensional case [17] and to study 
the collective modes in microbubble clouds [18] However, 
further work would be required to obtain quantitative 
predictions. 

Compressibility Tsunami waves generation by a mov- 
ing sea-bottom is on its own a difficult problem. Consid- 
ering the effect of compressibility in the generation of 
tsunami waves, the low resonant hydroacoustic mode of 
the water column at h = A/4 can affect the transmission 
of the sea bottom displacement to the free water surface 
via nonlinear mechanism depending on the moving bot- 
tom velocity. In addition, the hydroacoustics mode would 
constitutes the only tsunami generation mechanism in 
the absence of sea bottom residual displacement [19, 20]. 
A distribution of bubbles (cloud or metacloud) located 
at a significant depth [21] could affect this hydroacous- 
tic modes. However, for slow sea-bottom velocities and 



large generation areas compared with the water depth, 
the water compressibility can be neglected and then the 
generation of tsunami waves can be understood within 
the incompressible water framework[22, 23]. One may 
think that as we obtained our results from the acous- 
tic approximation and compressibility is a fundamental 
condition for acoustic waves to exist, metacloud and bub- 
bles in general have no chance to affect tsunami genera- 
tion when it is governed by incompressible mechanisms. 
However, due to a) the large differences between the com- 
pressibility of air and water and b) due to the deep sub- 
wavelength nature of the resonances studied, the incom- 
pressibility of water is also present in our proposal. The 
well known Rayleigh-Plesset equation, which is at the 
core of any study on air bubbles Moreover, in the spher- 
ical cloud model developed by d'Agostino and Brennen 
[12] incompressibility is recalled when deriving the scat- 
tering and absorption cross sections. We rederived their 
results in the linear regime including water compressibil- 
ity in the boundary conditions and the only difference 
between their derivation and ours lies on (i27rRi/\ + 1) 
factors. Provided that R\/\ « 1 (deep subwavelength 
regime) the incompressible approximation dominates. In 
simple words, as the wavelength is so large compared to 
the cloud (metacloud), the bubbles (clouds) only feel hy- 
drostatic pressure across the whole ensemble. Thus, in- 
compressible flow is not an issue that would turn off the 
collective metacloud resonance, although further research 
is certainly needed to know how the tsunami generation 
would be affected by these resonances. 

Nonlinearity It is well known that even for small am- 
plitude, bubbles can display a rich nonlinear behavior [6]. 
Cavitating clouds have also being studied mainly in the 
context of hydrodynamic systems [5, 6]. How nonlinear- 
ities would affect the collective oscillations of the meta- 
cloud and whether these nonlinearities would play or not 
a role in tsunami generation mitigation are challenging 
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questions that should be addressed in the future. 

Feasibility There are several factors which come into 
play concerning the feasibility of the metacloud approach 
for tsunami generation mitigation. If small bubbles were 
to be chosen, polydispersity should not be an issue since 
the lowest order modes are only slightly affected [18]. The 
generation of monodisperse microbubbles has been suc- 
cessfully achieved in laboratory environments [24], how- 
ever, another source of polydispersity will arise if the 
small bubbles are supposed to form clouds having tens 
of meters given by the static pressure difference along 
the cloud. Buoyancy and interaction of the clouds with 
ocean currents should also be considered. If instead of 
small bubbles, gas filled balloons ~ 1 m were to be cho- 
sen as the building block, buoyancy could be countered 
by means of ballast. In this case, the sphericity of the bal- 
loon would be compromised and its effect together with 
the effect of the covering membrane must be taken into 
account. Whether the cloud/metacloud should lie near to 
the sea surface or at a certain depth must be considered 
at the light of its influence on the tsunami generation. 
Although our proposal considers a certain degree of ran- 
domness within the metacloud, provided that a certain 
global filling fraction ci2 is reached and the clouds are 
at a minimum distance of each other (smaller than the 
interaction radius), in a real deployment the appropriate 
geometry might be an stochastic fractal one. 
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